figure(1);

lambda = 0.6;
theta = 0.1;

z = 2;
x = 0:0.01:3;

l = 0.5*(x - z).^2 + lambda*log(x/theta + 1);

plot(x, l);

% axis([min(x), max(x), min(l), max(l)]);